Termination of Multifractal behaviour for Critical Disordered Dirac Fermions 
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We consider Dirac fermions interacting with a disordered non-Abelian vector potential. The 
exact solution is obtained through a special type of conformal field theory including logarithmic 
correlators, without resorting to the replica or supersymmetry approaches. It is shown that the 
proper treatment of the conformal theory leads to a different multifractal scaling behaviour than 
initially expected. Moreover, the previous replica solution is found to be incorrect at the level of 
higher correlation functions. 
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There exists a lot of evidence that the wavefunctions 
in disordered systems exhibit multifractal behaviour near 
the localization-delocalization transition (see for exam- 
ple Multifractality manifests itself in the anoma- 
lous scaling behaviour of moments of the local density of 
states (LDOS) with system size. This anomalous scaling 
is described, for a d-dimensional system in a box of size 
L with short-distance cutoff a, by the spectrum r*(g) of 
exponents defined as 



(Mr)]"} ~ (a/Lr^ +d ; 



(1) 



Since r* (q) behaves linearly in q for a simple fractal (for 
example, a free wave in a box) , the nonlinear dependence 
of r*(q) signals the multifractality of the critical wave- 
functions. Such a behaviour has been observed in numer- 
ical calculations (for a review, see ||) and a parabolic 
dependence, r*(q) = (q — l)(d — aq), has been derived 
explicitly in a few different perturbative approaches such 
as the renormalization group (RG) one and the op- 
timal fluctuation- like methods ||. 

There are strong restrictions on the allowable be- 
haviour of T*(q) (see ||) which require among other 
things the termination of the parabolic dependence for 
large enough q (otherwise one ends up with unphysical 
arbitrarily negative exponents in (^)). One can imag- 
ine that this termination should be achieved when all 
multiple-loop corrections are taken into account in a per- 
turbative scheme. To check whether this hope is justified 
one inescapably needs to resort to nonperturbative meth- 
ods. 

In this respect two-dimensional disordered systems 
take a special place in investigating multifractality. On 
the one hand, there are critical disordered systems that 
are known to exhibit multifractality H|J: familiar ex- 
amples being quantum Hall systems and disordered sys- 
tems with spin-orbit interactions. On the other hand, 
in two dimensions one can study critical behaviour us- 
ing the powerful nonperturbative machinery of confor- 
mal field theory (CFT) 0. It is then natural to apply 
CFT to the problem of multifractality, and such attempts 



have indeed been initiated in Ref. by exploring ex- 
actly solvable models of Dirac fermions (introduced in 
JlOf in connection with the Integer Quantum Hall Effect) 
in Abelian or non-Abelian random vector fields. How- 
ever, the results obtained in these papers are somewhat 
puzzling. In particular the solution to the problem with 
SU(N) random vector field contains an infinite set of 
operators with negative conformal weights given by 
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leading via RG-type arguments to the set of (negative for 
large enough q) scaling exponents 



r*{q) = («-!) 2 



N-l 



(3) 



without explicit termination. The absence of the termi- 
nation mechanism should be taken more seriously here 
than in the case of RG calculations, since this CFT re- 
sult is nonperturbative and, if correct, exact. 

We reveal in this letter the existence of a termination 
mechanism for the non-Abelian case, emerging straight- 
forwardly from the correct full CFT treatment. This 
mechanism invalidates the result (0). For the sake of sim- 
plicity we illustrate this mechanism for the case N = 2, 
in which r* (g) becomes negative for the smallest value of 
q. Generalisation to arbitrary N is in principle possible. 

The termination mechanism we present is different 
from the one in the Abelian random vector potential 
case, which have been successfully calculated in Jll],[l2| by 
mapping into a random energy model or a type of Gaus- 
sian field theory. Moreover, for normalized wavefuntions, 
the termination of the multifractal spectrum can be ex- 
pected to originate from the instability of the effective 
Liouville field theory [§J for large q. 

The non-Abelian model has the additional virtue of be- 
ing exactly solvable without using replicas or supersym- 
metry (SUSY) Therefore we also get a chance to 
independently check the results of these two approaches 
P, ^4|p^| comparing them with our exact results. It is 



1 



well known that the replica approach fails to reproduce 
some (off-critical) nonperturbative results like the level 
statistics in random matrix theory |16|. For critical dis- 
ordered models one might expect more reliable results, 
since criticality implies universality. We here give an ex- 
ample demonstrating that this is not the case. Our exact 
solution shows that the higher correlation functions are 
identical with SUSY, but disagree with replicas. 

Another interesting aspect of our exact solution is in 
the nature of the CFT involved. One of the landmarks 
of conventional CFTs is the power-law dependence of the 
physical correlators. In fact, in many instances, CFTs are 
identified by just such a property. Recently, however, it 
has been discovered that logarithmic dependence in phys- 
ical correlators can appear in certain models outside of 
the usual class of so-called unitary minimal ones [fL7| . 
Logarithms are understood to be generated by degen- 
eracies in the spectrum of conformal dimensions of the 
theory: when two operators have dimensions becoming 
degenerate, they metamorphose into a logarithmic pair 
with unconventional correlators involving not only pow- 
ers, but logarithms. The model we are considering, in 
fact, enters this class of logarithmic CFTs. Moreover, as 
we shall see, it is the presence of such logarithmic oper- 
ators that provides us with the solution to the problems 
associated to negative scaling dimensions. 

Although the details of the calculation are somewhat 
involved, the conclusions can be relatively easily reached 
from our main result, which is the fusion rule (15). For 



the sake of clarity, though, we briefly outline in what 
follows the calculation process that leads to this result. 
Interested readers can find more extensive details in [jl8| . 

We consider N species of Dirac fermions living in a 
2+1-dimensional space and interacting through a dis- 
ordered vector potential transforming like the ad- 
joint of an su(N) algebra A, to which they are cou- 
pled minimally. The disorder allows for hopping between 
the different species. Since the vector potential is time- 
independent, different Matsubara frequencies do not cou- 
ple, and can be treated independently by a Euclidean 
two-dimensional theory with explicit frequency depen- 
dence. In fact, for a given realization of the disorder, the 
partition function takes the form of the fcrmionic path 
integral with the Dirac action 



S[*,a;,A„] = / d 2 x^(x)[I 



iui + iA]^ (4) 



(since we are in a two-dimensional Euclidean space, we 
take the Pauli matrices as Dirac 7 matrices). 

Arbitrary products of disorder-dependent single- 
particle Green's functions can then be calculated and 
averaged over the vector field distribution functional 



representing the usual 5-correlated Gaussian white noise 
for the random vector potential. 

In the limits of infinite disorder strength g — > 00 and 
of vanishing frequency u> — > 0, the theory becomes con- 
formally invariant. Correlators can then be calculated 
according to the general principles behind CFT. 

The derivation involves the following key moments. 
We separate the fermionic action into chiral parts, using 
holomorphic and antiholomorphic derivatives and fields 
(2d = cL,2<9 = d+), A± = Ai ± iA 2 . It is important 
to note that now A± 6 su c (N), the complex exten- 
sion A c of A. We then parametrize the vector fields 
by fields g± belonging to the complex extension G c of 
the group G = SU{N) as A±(x) = id±g±(x)g± 1 (x). 
The reality condition A\(x) = A-(x) translates into 
g\ (%) = 9-( x )- We will use the notation g+(x) — g(x). 
This reparametrization induces a non trivial Jacobian in 
the path integral @Jl§] 



VA. = VA1VA2 = VG c e 2NW ^^ 



(6) 



where VG C is the Haar measure over SU C (N), and 
Wfgtg] is the Wess-Zumino-Novikov-Witten functional 
here on level k = —27V, for the field combination 
_ 'g. This is a well-known functional for which correla- 
tors are in principle known po[ . 

The next step is to decouple the fermions from the 
random potential by the transformations 

* ±(x) - g±(x)9' ± (x), *l(x) - ^(x)g^ x (x) (7) 

The Jacobian for this decoupling has the very important 
property of being proportional to the partition function 
at fixed disorder |]l9|| , thus cancelling it when computing 
the correlations for a given realization. This removes the 
need to invoke either the replica or SUSY methods to 
perform explicitly the disorder averaging. 

Let us then consider correlators of the local operator 



M(z, z) = Tr** = *'Jh ab *' + 9' + l h **- 



(8) 



P[A„] = I J (fxTrA^A^x) 



(5) 



(in which h a i, = [g^g] a b, with th SU(N) indices) which 
relates to the LDOS of the system and couples to the fre- 
quency perturbation in the original action. Notice that 
by doing this, we explicitly sum over the SU (N) indices: 
the product g^ g is invariant under left-multiplication of 
g by u G SU(N). This is a crucial point: g^g does not 
live on the complexified group manifold of SU (N), but 
rather on the coset space SU C (N)/SU{N). The SU(N) 
path-integration then simply factorizes out of the effec- 
tive generating functional. 

In general, when one deals with CFTs, physical states 
are associated to states in a highest-weight representa- 
tion of the Virasoro algebra 0. Quantum mechanical 
commutation rules are replaced in the radial quantiza- 
tion formalism by so-called fusion rules, which state the 
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short-distance singular behaviour of products of opera- 
tors in the complex plane defined by the Euclidean space- 
time variables z = x\ + ix 2 , z = x\ — ix 2 . Such fusion 
rules are in principle completely obtainable from the four- 
point correlation functions of the physical operators, for 
which very powerful and extensive calculational methods 
are known. 

One of the crucial aspects of the derivation is the proof 
that the WZNW model of level k = -2N on the coset 
space SU (N)/SU(N) actually carries a representation 
of the Virasoro algebra for SU(N) for the same (analyti- 
cally continued to negative) level. We refer the reader to 
fl8f for an explicit proof of this statement. 

The main procedures for obtaining the four-point coset 
correlators 



with 

H tJ (z,z) = a\z(l - z)f 2 [F?{z)F*(z) + (a <-> b)] 



H 



(h aibl (z 1 ,z 1 )h b2 \ 2 (z 2 , z 2 )h a3b3 {z 3 , Z 3 )h b ^{zi, Zi)) (9) 



and thus the fusion rules, in the case of WZNW mod- 
els is very clearly set out in [^0|. The "conformal 
bootstrap" provides differential equations (the Knizhnik- 
Zamolodchikov equations) for the so-called conformal 
blocks, whose solutions turn out to be hypergeometric 
functions of the variable z. The peculiarity that we en- 
counter in the ensuing CFT, i.e. the one with SU(2) 
Virasoro algebra on level k = —4, is that some conformal 
blocks have logarithmic behaviour. In fact, these confor- 
mal blocks read 

F 1 a (z) = (l-z)F(3/2,5/2;2;z) 



Ft{z) = 



4/3 



+ F?(z) Inz- 4/3 + (1 - z)K n (z) 



F2(z)=zF(3/2 1 5/2;3;z) 



F b 2 {z) 



1 

2^ 



16/3 



+ F2(z)lnz-4/3 + zK 12 (z) 



(10) 



where F(a, 6; c; z) is the hypergeometric function, and 
K\\,K\ 2 are some functions regular as z — > 0. 

Enforcing the locality conditions of quantum field the- 
ory then translates into "gluing" together the confor- 
mal blocks for the spacetime variables z and z in such 
a way that the combination is single-valued in the com- 
plex plane defined by the variable z (we refer the reader 
to plj] for the basic explanations). The main point that 
we want to stress here is that these gluing conditions 
in logarithmic CFTs are not, as was once thought, un- 
solvable, but rather a bit peculiar (for more details, see 
fyl). The crucial observation to make is that, as we take 
z — > e 2m z, only the combination hi zz remains invariant. 
This dictates the proper gluing procedure. 

Leaving out the explicit details, we here simply pro- 
vide the final expression for the coset correlator (z — 
Z12Z34 . — 
Z13Z24 ' y 



Zj): 



I2 — & €110,460,20,3 (12) 



and for the resulting fusion rule between the coset oper- 
ators (a similar one has been obtained also in Ref. p2|): 



r TTh(zi,Z 1 )h- 1 (z 2 ,Z 2 ) ~ M 3 / 2 



\[^A{z 2 )+ 

■i Zi 2 



—A(z 2 )} + [41 + 2D(z 2 , z 2 ) + In \z 12 \C{z 2 , z 2 )\ 
z\ 2 



(13) 



in which X is the identity operator. The correlators of 
the operators appearing in (O) are 



H = \z 13 z 24 



3/2 



E 



IiIjHij(z,z) 



(A( Zl )A(z 2 )) ~ z( 2 ; (A( Zl )A(z 2 )} ~ z( 2 
{D(z uZl )D(z 2 ,z 2 )) - -ci -ln|zi 2 | 
(D{z 1 ,z 1 )C{z 2 ,z 2 )) ~ 1; 

(C(z 1 ,z 1 )C(z 2 ,z 2 )} = (14) 

where ci are some constants, unimportant for our pur- 
poses. 

There are many remarks that we can make from 
this unconventional Operator Product Expansion (OPE). 
First and foremost, wee notice that the most relevant op- 
erators appearing in the OPE (|l^), A(z) and A(z), pos- 
sess conformal weights (—1,0) and (0,-1) respectively. 
Usually, the fusion rules for WZNW models [^0| would 
imply that the adjoint operator, whose conformal weights 
are (—1,-1), should appear in the OPE (13). But the 
term in the four-point function pointing to the presence 
of such an operator, does not appear since we are not 
allowed, by the requirement of single valuedness, to mul- 
tiply the logarithmic solutions in the holomorphic and 
antiholomorphic sectors together. This requirement can- 
not be seen from the chiral conformal algebra studied 
in Q: it can only come out of the solution for the full 
correlator that we have obtained. 

The four-point M. correlator can be calculated along 
the same lines as the coset correlator ( |iT| ) by including 
the free fermion contributions and contracting the indices 
appropriately. This four-point function allows us to ex- 
tract the OPE of M. with itself, which in turn determines 
the scaling of the local moments. We find that the correct 
OPE reads M: 



M{l)M{2) 



1 



\Z\2 



1/2 



J + £(2) + iln|z 12 |C(2) 



(15) 



(11) where a is some constant. The crucial fact is that the op- 
erator in the symmetric representation does not appear 
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here again, like in (|13|). D and C operators fuse as (14) 
whereas M. , D and C fuse as 



D{l)M{2) 
C(1)M(2) 



• [In |*i2 1 
M{2) 4 



c 2 ]M(2) 



(16) 



The fusion rules (|15|,[16|) constitute the most impor- 
tant result of our paper. The remarkable feature of these 
OPEs is the appearance of logarithmic operators D and 
C. Their origin lies in logarithmic singularities present 
in multi-point correlation functions. These operators do 
not form the usual diagonalizable representations of the 
Virasoro energy operator Lq |17| which is evident from 
their unconventional correlation functions ( |l4| ) . 

The first comment that we can make about the OPE 
( |j~5| ) is that it invalidates the previous replica method 

FT 



treatment |14 
M{l)M{2) 



for which the corresponding OPE read 
1 



\Z\2 



1/2 



[1 + z[D{2) +hx\z\C{2)]- 
±z{D(2)+ln\z\C[2)] + .. 



(17) 



in which C, D, C, D were chiral logarithmic pairs with 
conformal dimensions (1,0) and (0,1) respectively. This 
is different from the pairs introduced in ( |l4j ) which have 
dimensions (0,0). Even though the two-point functions 
coincide in both treatments, the higher-point functions 
are different. Operator dimensions are calculated cor- 
rectly by replicas, but correlators of higher order are not. 

Our result ([if]) will also lead to a markedly different 
behaviour than the one associated to equation (||) . Since 
the LDOS p is related to the imaginary part of A4, its lo- 
cal moments p q will scale at most like the most relevant 
part of A4 q , which can be obtained by point-splitting 
from ([15]). We can draw the following comparison table 
between our results for the conformal weights h q of the 
most relevant operator contained in A^' as a function 
of the power q, and the ones obtained in the previous 
treatment [HI: 



q 


1 


2 


3 


4 




h q ( [9]) 


1/8 


- 1 


-9/4 


-4 




h q (exact) 


1/8 





1/8 








The parabolic series of negative exponents is seen to be 
cut right from the beginning for SU{2). Note that log- 
arithms are weight-zero objects, so in fact they appear 
implicitly in our weights above for all q. The presence of 
these logarithmic prefactors makes the usual treatment 
with r* exponents somewhat unsatisfactory, since it can 
only describe pure power scaling and not the logarithmic 
corrections to it that we have found. 

In conclusion, we have shown that moments of the lo- 
cal density operators in the problem of non-Abelian ran- 
domness obey different scaling relations than expected 
P|, that do not fit in the standard multifractal descrip- 
tion. The negative-dimensional operators are suppressed 



by the presence of logarithmic operators in the relevant 
OPEs. Moreover, although replicas give correct primary 
field scaling dimensions, they fail to reproduce the de- 
tailed form of the correlators. Further details can be 
found in [[§]. 
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